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Abstract: Approximately taking into account the higher-order effects due to multi- 
ple-gluon initial-state radiation, we extract from the latest Tevatron data of prompt 
J/ip hadroproduction the leading colour-octet matrix elements within the nonrelati- 
vistic-QCD (NRQCD) factorization formalism. We find that the matrix elements 
which describe the formation of J/if> mesons from colour-octet cc pairs in the 
*0 . angular-momentum states 2S+l Lj = and 3 Pj, with J = 0, 1, 2, which are respon- 

^ sible for the excess of the predicted cross section of inelastic J /if) photoproduction 

\ over the existing HERA data at high values of the inelasticity variable z, are sig- 
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nificantly reduced. As a consequence, the Tevatron and HERA measurements of 



Qs ' inclusive J /if) production are reconciled in the NRQCD framework. 

a: 

P> 1 Introduction 

The cross section of inclusive J / if) hadroproduction measured in pp collisions at the Fermilab 
Tevatron |2|] turned out to be more than one order of magnitude in excess of what used to 
be the best theoretical prediction, based on the colour- singlet model (CSM). As a solution to 
this puzzle, Bodwin, Braaten, and Lepage || proposed the existence of so-called colour-octet 
processes to fill the gap. The central idea is that cc pairs are produced at short distances 
in colour-octet states and subsequently evolve into physical (colour- singlet) charmonia by the 
nonperturbative emission of soft gluons. The underlying theoretical framework is provided 
by nonrelativistic QCD (NRQCD) endowed with a particular factorization hypothesis, which 
implies a separation of short- distance coefficients, which are amenable to perturbative QCD, 
from long-distance matrix elements, which must be extracted from experiment. This formalism 
involves a double expansion in the strong coupling constant a s and the relative velocity v of 
the bound charm quarks, and takes the complete structure of the charmonium Fock space into 
account. 

In the case of inelastic J /if) photoproduction, NRQCD with colour-octet matrix elements 
tuned || to fit the Tevatron data H predicts || at leading order (LO) a distinct rise in cross 
section as z — > 1, where z is the fraction of the photon energy transferred to the J /if) meson 
in the proton rest frame, which is not observed by the HI [B] and ZEUS |7| collaborations at 



1 To appear in the proceedings of the workshop on Monte Carlo Generators for HERA Physics, DESY, 
Hamburg, 21 April 1998 - 5 February 1999, edited by T. Doyle, G. Grindhammer, G. Ingelman, and H. Jung. 



1 



HERA. This colour-octet charmonium anomaly has cast doubts on the validity of the NRQCD 
factorization hypothesis 0, which seems so indispensible to interpret the Tevatron data in a 
meaningful way. 

Here, we report on an attempt || to rescue the NRQCD approach by approximately taking 
into account dominant higher-order (HO) QCD effects. The basic idea is as follows. The 
predicted excess over the HERA data at z close to unity is chiefly generated by colour-octet 
cc pairs in the states ^o, 3 -Po ; an d 3 P2 @, where we use the spectroscopic notation 2S+1 Lj to 
indicate the spin S, the orbital angular momentum L, and the total angular momentum J. On 
the other hand, in hadropro duct ion at the Tevatron, the contributions from the colour-octet 

and 3 Pj states fall off much more strongly with increasing transverse momentum (p?) than 
the one due to the colour-octet 3 Si state j|, which is greatly suppressed in the quasi-elastic 
limit of photoproduction [Q. Consequently, the nonperturbative matrix elements which are 
responsible for the colour-octet charmonium crisis are essentially fixed by the Tevatron data in 
the low-pr regime. This is precisely where the LO approximation used in Ref. @] is expected to 
become unreliable due to multiple-gluon radiation from the initial and final states. In Ref. |9|], 
this phenomenon was carefully analyzed in a Monte Carlo framework and found to significantly 
increase the LO cross section. In Ref. ||, fits to the latest prompt J/ip data taken by the 
CDF collaboration at the Tevatron were performed incorporating this information |J on the 
dominant HO QCD effects. The resulting HO-improved NRQCD predictions for inelastic J/i/j 
photoproduction at HERA do not overshoot the HI |J and ZEUS @ data any more. 



2 Theoretical input 

The underlying theoretical framework is explained in Ref. ||. If is of order Mj/^ or below, 
we adopt the fusion picture, where the cc bound state is formed within the primary hard- 
scattering process. In the high-pr regime, we work in the fragmentation picture, where the cc 
bound state is created from a single high-energy gluon, charm quark or antiquark which is close 
to its mass shell. We take the renormalization scale \i and the common factorization scale Mf to 



be jj, = Mf = rriTi where my = y4m^ + p\ is the J /if) transverse mass. We define the starting 
scale no of the fragmentation functions (FF's) as fio — 2m c = Mj/$. For our LO analysis, 
we choose CTEQ4L [TIJ and GRV-LO |TT| as the proton and photon PDF's, respectively, and 



evaluate a s from the one- loop formula with = 236 MeV |IJ . Whenever we include higher 



orders, we adopt the MS renormalization and factorization scheme and employ CTEQ4M |10 



GRV-HO 0, and the two-loop formula for a s with = 296 MeV |TJ. 



MS 

Unfortunately, not all ingredients which would be necessary for a fully consistent NLO anal- 
ysis are yet available. In the case of fusion, the NLO corrections to the partonic cross sections 



are only known for direct photoproduction in the CSM Furthermore, in the case of frag- 
mentation, the NLO corrections to the FF's at the initial scale fio are still unknown. In the 
case of direct J / if) photoproduction under typical HERA conditions, the QCD correction factor 



K to the inclusive cross section in the CSM was found |12| to be as low as 1.2 in the inelastic 
regime z < 0.9. It is plausible that the K factors for the color-octet and resolved-photon pro- 
cesses should be modest, too. However, the situation should be very different for inclusive J/ip 
hadroproduction at the Tevatron, especially in the low-pr range, where one expects substantial 
HO QCD effects due to multiple-gluon radiation. Such effects were estimated for the fusion 
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Table 1: Values of the J/ip matrix elements resulting from the LO and HO-improved fits to 
the CDF data [§. Here, = (C J/ ^[8, ^o]) + (r/m 2 c ) (c J ^[8, 3 P ]). The first 10 (last 2) 

of the 11 data points are described in the fusion (fragmentation) picture. 





LO 


HO 


{0 J /1'[l, a S 1 ]\ 

r 
Xdf 


(7.63 ± 0.54) • lO" 1 GeV 3 
(3.94 ±0.63) ■ lO- 3 GeV 3 
(6.52 ± 0.67) • 10~ 2 GeV 3 
3.47 
7.49/12 


(1.30 ±0.09) GeV 3 
(2.73 ±0.45) • 10- 3 GeV 3 
(5.72 ± 1.84) • 10~ 3 GeV 3 
3.54 
3.96/12 



mechanism in Ref. |9| by means of the Monte Carlo event generator PYTHIA [T3J after im- 
plementing therein the relevant colour-octet processes, and they were indeed found to be very 
sizeable. The impact of these effects on the fit to the latest CDF data is demonstrated in 
Table § 



3 Predictions for charmonium photoproduction 

We now explore the phenomenological consequences of this HO improvement for inclusive J/i/j 
photoproduction in ep collisions at HERA, with beam energies E e = 27.5 GeV and E p = 
820 GeV in the laboratory frame, assuming the maximum photon virtuality to be Q max = 

4 GeV 2 . As in the HI |6j and ZEUS 0] publications, we convert the ep cross sections to averaged 
7P cross sections by dividing out the photon-flux factor. The contribution due to ip' mesons 
with subsequent decay into J/ip mesons is approximately taken into account by multiplying 
the theoretical predictions by an overall factor of 1.15. The data are mostly concentrated in 
the low-pr range, where the fusion picture should be valid. 

In Fig. 1, we compare our LO and HO-improved predictions for the p^, z, and y distri- 
butions with the HI || and ZEUS Qj data. Here, y is the J/ip rapidity in the laboratory 
frame, which is taken to be positive in the proton flight direction. The circumstance that 
(O j ^[8, %}) and (C J ^'[8, 3 P ]) are not separately fixed by the fit to the CDF data § induces 
some uncertainty in the colour-octet contributions to the cross sections of direct and resolved 
photoproduction and thus also in the total cross section. This uncertainty is encompassed by 
the results for ((9 J ^[8, %]) = M J r ^ and (0 J ^[8, 3 P ]) = and those for (C J ^'[8, %)) = 

and (^O j ^[8 } 3 P ]) = (m 2 /r)M r J// ^, which are actually shown in Fig. 1. We observe that, at LO, 
the colour-octet contribution of direct photoproduction is dominant for z > 0.5. Thus, it also 
makes up the bulk of the p\ and y distributions, which are integrated over 0.4 < z < 0.9. This 
contribution is responsible for the significant excess of the LO predictions over the experimen- 
tal results for dojdp\ and do jdz at low px and high z, respectively. On the other hand, the 
HO-improved predictions tend to undershoot the data leaving room for a substantial K factor 
due to the missing NLO corrections to the partonic cross sections. Now, the colour- singlet 
contribution of direct photoproduction, which is well under theoretical control [12|, is by far 



3 



dominant, except in the corners of phase space, at z < 0.15 and z > 0.85, where the colour- 
octet contributions of resolved and direct photoproduction, respectively, take over. Of course, 
we should also bear in mind that the predictions shown in Fig. 1 still suffer from considerable 
theoretical uncertainties related to the choice of the scales /i and Mf, the PDF's, and other 
input parameters such as m c and A^ 4 ^ ||12|| . From these observations, we conclude that it is 
premature at this point to speak about a discrepancy between the Tevatron || and HERA 
H |7| data of inclusive J/ip production within the framework of NRQCD [HI. 



4 Conclusions 

We determined the J/ip colour-octet matrix elements which appear in the NRQCD expansion 
[|3| at leading order in v by fitting the latest Tevatron data of prompt J/ip hadroproduction 0. 
We found that the result for the linear combination M; 7 ^ of ((9 J/ ^[8, ^o]) and (0 W [8, 3 P O [) 
is substantially reduced if the HO QCD effects due to the multiple emission of gluons, which 
had been estimated by Monte Carlo techniques |J, are taken into account. As an important 
consequence, the intriguing excess of the LO NRQCD prediction for inelastic J/ip photopro- 
duction at z close to unity |J over the HERA measurements [|, [7| disappears. We assess this 
finding as an indication that it is premature to proclaim an experimental falsification of the 
NRQCD framework on the basis of the HERA data. Although we believe that our analysis 
captures the main trend of the HO improvement, we stress that it is still at an exploratory level, 
since a number of ingredients which would be necessary for a fully consistent NLO treatment 
of inclusive J/ip hadroproduction and photoproduction are still missing. 
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Figure 1: The HI ^ and ZEUS /1W data on inelastic J /if) photoproduction are compared with the 
LO and HO-improved predictions. The total results (solid lines) are built up by the direct-photon 
colour-singlet (dotted lines), direct-photon colour-octet (dashed lines), resolved-photon colour- 
singlet (dot-dashed lines), and resolved-photon colour-octet (dot- dot- dashed lines) contributions. 
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